https://doi.Org/10.38032/jea.2020.02.001 


4 SciEn 

Journal of Engineering Advancements Vol. 01(02) 2020, pp 29-36 


Numerical Analysis of the Aerodynamic Characteristics of NACA-4312 Airfoil 

Md Rhyhanul Islam Pranto*, Mohammad Ilias Inam 
Department of Mechanical Engineering, Khulna University of Engineering & Technology, Khulna-9203, BANGLADESH 
Received: May 12, 2020, Revised: June 17, 2020, Accepted: June 18, 2020, Available Online: June 24, 2020 


ABSTRACT 


The aim of the work is to investigate the aerodynamic characteristics such as lift coefficient, drag coefficient, pressure distribution over 
a surface of an airfoil of NACA-4312. A commercial software ANSYS Eluent was used for these numerical simulations to calculate 
the aerodynamic characteristics of 2-D NACA-4312 airfoil at different angles of attack (a) at fixed Reynolds number (Re), equal to 
5x10^. These simulations were solved using two different turbulence models, one was the Standard k — s model with enhanced wall 
treatment and other was the SST k — o) model. Numerical results demonstrate that both models can produce similar results with little 
deviations. It was observed that both lift and drag coefficient increase at higher angles of attack, however lift coefficient starts to reduce 
at a =13° which is known as stalling condition. Numerical results also show that flow separations start at rare edge when the angle of 
attack is higher than 13° due to the reduction of lift coefficient. 
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1. Introduction 

Aerodynamics is defined as a study of the resulting effects of 
relative motion between air molecules and body surfaces 11]. 
Understanding the air flow behavior around a moving body is 
very important because we can control the air flow to our 
advantage in the case of aircraft, wind turbine, drone, and many 
more. Aerodynamics focuses on studying this phenomenon by 
applying basic physics laws such as Newton's laws and the 
Navier-Stokes equation. To study the aerodynamic profile of a 
particular object, it is necessary to define the object's shape, 
where the airfoil or the wing cross-section is used in the case of 
an airplane. The cross-section of the typical aircraft wing is an 
airfoil and is largely responsible for producing the 
forces that sustain the aircraft in flight 12]. Airfoil is known as 
the cross-sectional shape of a wing, blade (of a propeller, rotor, 
or turbine) which is placed in an airstream in order to generate 
useful aerodynamic forces. An airfoil's aerodynamic profile can 
be examined in two ways, one experimentally and other 
computationally. The aerodynamic characteristics of an airfoil 
are obtained in the experimental method by placing the airfoil (or 
wing model) in a controlled wind tunnel and recording the 
velocity and pressure distribution around the airfoil. The same 
experiment can be done in the computational approach using a 
CED code. Both approaches are usually complementary to each 
other and can help us design a specific purpose for the airfoil 13]. 

In the early 1800s, it was discovered that a curved surface 
produces more lift than a similar size fiat plate by Sir George 
Cayley 14]. The most effective way to do this is to use an airfoil. 
In the Langley two-dimensional low turbulence tunnels tests a 
considerable amount of airfoil data has been accumulated by 
National Advisory Committee for Aeronautics (NACA) in the 
United States. The development of NACA airfoils which are now 
in common use was started in 1929 with a systematic 
investigation of a family of airfoils in the Langley variable- 
density tunnel 15]. Airfoils of this family are designated by 
numbers having four digits. Experiments were performed to test 


different airfoils and to determine the lift and drag coefficient for 
different airfoils and to report their results. Jacobs et al. (1933) 
performed a wide variety of experiments comparing the 
geometric airfoil parameters and angle of attack and wind speed 
for different NACA airfoils and published the results of the lift 
and drag coefficient for 78 different NACA airfoils 16]. Since 
these experiments were carried out in a wind tunnel with a 
relatively high free-stream turbulence intensity of about 2%, the 
results are expected to vary from low free-stream turbulence 
conditions. NACA airfoils have been used extensively as 
aerodynamic test models. NACA-0012 and NACA-4412 
profiles are arguably the most studied airfoils. Ravi et al. 17] 
predicted numerically a transition model of an incompressible 
laminar to turbulent flow over NACA-4412 airfoil at Reynolds 
number of 3 X 10^. Eleni et al. 18] also carried out studies on 
the variation of lift and drag coefficients from flow around 
NACA-0012 airfoil at Reynolds number of 3 X 10^ for different 
turbulence models. In addition to the results available for a 
NACA-0012 and NACA-4412 airfoil profile, a substantial 
experimental or numerical database is not available for NACA- 
4312 airfoil. This study is motivated by the need to complement 
the presently limited body of knowledge for NACA-4312 airfoil. 
Thus, the present investigation is focused on examining the effect 
of the Reynolds number and the angle of attack on the 
performance characteristics of a NACA-4312 airfoil and relating 
the performance characteristics. 

In this study, CED approach is used to determine the 
aerodynamic profiles of the NACA-4312 airfoil. The CED 
packages contain three main elements which are Pre-processor, 
Solver and Post-processor 19]. Recently CED has been the 
method of choice in the aerospace, automotive and many 
industrial components. CED is vastly used in the field of 
aerodynamics because it is cheaper than experimental process 
and also gives more accurate results. The transition from the 
laminar to the turbulent flow phase plays a very important role in 
simulating the flow over an airfoil. The fluid flow over the airfoil 
exerts a pressure force perpendicular to the upper and lower 
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surfaces along with shear force which is parallel to the surfaces. 
The resultant of these two forces is the area of interest. The 
normal component of the force is known as lift force and the 
force which is acting in the flow direction is known as drag [10]. 

It was observed from the previous literature review that there 
is a lack of understanding of the characteristics of NACA-4312 
airfoil. In this paper a series of numerical simulations were 
carried out to analyze the characteristics of NACA-4312 airfoil 
with cord length of 1000 mm and Reynolds number of 5 X 10^. 
ANSYS Fluent was used to solve the steady-state RANS 
(Reynolds Average Navier-Stokes) equation with two different 
turbulence models. 

2. Methodology 

2.1 Airfoil 

In this paper, NACA-4312 airfoil was selected for the 
simulations, the profile is shown in Fig. 1. The NACA-4312 
airfoil means that it has a maximum camber equal to 4% of the 
cord which is located at 30% of the cord from the leading edge 
with a maximum thickness of 12% of the chord. In this paper, 
cord length was considered 1000 mm. Airfoil was created by 
ANSYS design modeler by importing coordinate file of 
NACA-4312. 
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Fig. 1 Geometry of NACA-4312 airfoil. 

2.2 Computational Method 

These simulations were conducted by ANSYS Fluent. The 
problem was solved in steady-state with two turbulence 
models. These simulations were conducted at fixed Reynolds 
number (Re), equal to 5x10^, to show the effect at the 
transition region [11]. Air was assumed as working medium 
with a constant density (p) and viscosity (p), whereas p = 
1.225 kg/m^, and p = 1.7894 x 10“^ kg/ms [12]. The steady- 
state RANS (Reynolds Average Navier-Stokes) equation was 
solved using the Green-Gauss cell based gradient option and 
pressure based solver was selected as the flow is 
incompressible. The RANS equations are time- 
averaged equations of motion for fluid flow, as bellows: 
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where p indicates density, u indicates velocity and ju 
indicates dynamic viscosity of the fluid. The left side of this 
equation describes the fluid element's mean momentum change 
to the flow instability and convection by the mean flow. This 
adjustment is regulated by the mean body force, the isotropic 
stress due to the mean pressure field, the viscous stresses, and 
visible stress due to the fluctuating velocity field, commonly 
referred to as the Reynolds stress (—pzrjzrj). This nonlinear 
stress term requires additional modeling to solve the RANS 
equations and has resulted in many different turbulence models 
[13]. 

In this paper, the Standard k — s model and the SST k — 
0 ) model are used to predict the effects of turbulence in flow 
over the airfoil. 

2.2.1 Standard k-s Model 

K-epsilon (k — s) turbulence model is the most common 
model used in Computational Fluid Dynamics (CFD) for 
simulating flow characteristics for turbulent flow conditions. It 
is a two-equation model, which provides a general definition of 
turbulence via solving two transport equations (PDFs) one for 
the turbulence kinetic energy (k) and the other is its dissipation 
rate (e) [13]. The turbulent kinetic energy k and dissipation £ 
are obtained from the following transport equations (3) and (4): 
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In these equations, 

Gk= Generation of turbulence kinetic energy due to mean 
velocity gradient 

YM=Fluctuation in compressible turbulence to overall 

dissipation rate 

Cis> ^38 = Constants 

(Jk> ^ 8 = Turbulent Prandtl number for k and £ 

Sk,S^ = User defined source terms 

Here turbulent viscosity = pC^ — ; where is C^a 
constant. 

2.2.2 SST k-co Model 

In CFD analysis SST (Shear Stress Transport) k — o) is sl 
widely used and two-equation eddy-viscosity turbulence 
model. It provides a general solution for the turbulence kinetic 
energy (k) and specific dissipation rate of eddy viscosity (o)). 
This model gives more accurate results in numerical flow 


30 













M.R.I Pranto and M.I. Inam /JEA VoL 01(02) 2020, pp 29-36 


analysis of airfoil than the standard k — o) turbulence model as 
it includes transportation of turbulence sheer stress. The proper 
transport behavior can be obtained to the formulation of eddy- 
viscosity [14] .The turbulent kinetic energy k and dissipation 
rate of eddy viscosity o) are obtained from the following 
transport equations (5) and (6): 

(5) 
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Gk= Generation of turbulence kinetic energy due to mean 
velocity gradient 
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where S is the 


strain rate magnitude, F 2 is a constant, and 
oc=Angle of attack. 


2.3 Boundary Condition with domain 

For these simulations, a computation domain was created 
around the NACA-4312 airfoil, shown in Fig. 1. Chord length 
of NACA-4312 airfoil was assumed to be equal to 1000 mm. 
To minimize the boundary effect, the domain was extended 
12.5C (C indicated chord length) along the upstream and 20C 
along downstream from the trailing edge, like Fig. 2 [15]. In 
the domain, airfoil surface was assumed no-slip boundary 
condition, BAFED was assumed constant velocity inlet and 
pressure outlet at BCD. For this fixed Re, the inlet velocity u, 
was assumed to 7.5 m/s. In these simulation angle of attack (a) 
was changed by changing flow direction instead of rotating 
airfoil which produced the same effect on the airfoil. The 
velocity components along x and y direction were calculated by 
u cos a and u sin a, respectively, for different angle of attack 
(a). 


2.4 Mesh generation and Wall Treatment 

The C-type structural mesh was created for better 
convergence and control of the wall function, shown in Fig. 3 
and Fig. 4. Fig. 3 depicts the structural meshing with 
quadrilateral elements which was done by ANSYS Meshing 
whereas Fig. 4 illustrates the meshing quality around the airfoil 
profile. 

The application of wall functions near the wall region may 
significantly reduce the processing and storage requirements 
while producing an acceptable degree of solution. To ensure 
sufficient boundary layer modeling the inflation was set to 1.15. 

The effectiveness of k — £ model is beyond = 30. The 
non-dimensional wall parameter is defined as: 


A Velocity Inlet B 




Fig. 3 Mesh among the whole domain. 



Fig. 4 Mesh around the NACA-4312 airfoil. 
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+ yJ(TCD/ p) 

y =yx - 

M 

Here is the distance from the wall to the centroid of the 
first fluid cell. It was observed for all simulations that the value 
of y'^was higher than 30 and lower than 60. 

2.5 Mesh Independence Test 

A series of simulations were conducted for mesh 
independency tests. A different number of meshes were created 
by diving circular and rectangular sections with different 
numbers. Fig. 5 depicts the variation of lift coefficient for 
different numbers of elements of the mesh for an angle of attack 
5°. It is observed from figure that meshes with higher than 
105000 number of elements can produce accurate results with 
minimum deviation. As a result mesh with element number 
105000 were selected for further simulation. 



Fig. 5 Variation of coefficient of lift with no. of element. 

3. Results and Discussion 

Fig. 6 depicts the variation of lift coefficient for different 
angles of attack for two different turbulence models. It is 
observed that the lift coefficients are almost equal for both 
models up to a = 6®, however, the SST k — co model predicts 
higher lift coefficient compared to the standard k — s model 
when a > 6®. Fig. 6 also shows that the lift coefficient increases 
with angle of attack, for both models, up to 12®. Later on, it 
decreases with a due to the flow separation at the trailing edge, 
which indicated stalling occurs at the range of 12®13®. From 
0° angle of attack to 12° angle of attack the lift curve is almost 
linear. 

Fig. 7 shows the effect of the angle of attack (a) on the 
drag coefficient for two different turbulence method. It is 
observed that the drag coefficient estimates higher values with 
the standard k — s model compared to the S ST k — o) because 
of the transportation of the shear viscosity. For both models, the 
drag coefficient increase with a, however, the increment rate is 
higher at higher a due to the flow separation at trailing edge. 
The drag coefficient still increases after the stalling angle 
however lift coefficient decreases. 



Fig. 6 Lift coefficient vs angle of attack. 



Fig. 7 Drag coefficient vs angle of attack. 



Fig. 8 Pressure coefficient on the airfoil surface at 10° angle 
of attack. 
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Fig. 9 Contours of pressure with the Standard k-s model for different angles of attack of NACA-4312 airfoil. 
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Fig. 10 Contours of pressure with the SST k-co model for different angles of attack of NACA-4312 airfoil. 
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Fig. 11 Velocity contours with the Standard k-s model 
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Fig. 12 Velocity contours with the SST k-co model 


34 




































M.R.I Pranto and M.I. Inam /JEA VoL 01(02) 2020, pp 29-36 


Fig. 8 demonstrates pressure distribution on the airfoil 
surface for both models at an angle of attack a = 10®. The upper 
and lower portion of the graph represent pressure coefficient 
distribution on the lower surface and upper surface of the 
airfoil, respectively. It is observed that the upper surface have 
lower pressure and lower surface have higher pressure which 
produces lift. Maximum lower and higher pressure observed 
near to the leading of the lower and upper surface, respectively, 
which starts to reduce along the trailing edge. 

3.1 Pressure Contours 

Fig. 9 and Fig. 10 presents pressure contours for the 
Standard k — s turbulence model and the SST k — co 
turbulence model, respectively, at different angle of attack. In 
these two figures, the effect of angles of attack is shown by 
changing the flow directions instead of rotating airfoil, which 
gives the same effect on airfoil. Similar pressure distribution 
pattern is observed for the Standard k — e and SST k — o) 
turbulence model, present in Fig. 9 and Fig. 10, respectively at 
lower angle of attack, whereas slightly varies in magnitude at 
higher angle of attack though the pattern of pressure 
distribution is similar. It is clearly observed from both figures 
that the upper surface has lower pressure and the lower surface 
has higher pressure. It is also observed from both figures that 
negative pressure region on the upper surface which is situated 
almost around the whole surface which starts to decrease with 
the increase of angle of attack, however pressure on the lower 
surfaces increases with angle of attack. Though pressure 
difference between lower and upper surface increase until flow 
separates on the upper surface at a certain angle of attack. The 
lift coefficient increases in both turbulence model with the 
increase of angle of attack and it reaches the stalling position 
near 13° angle of attack. Between the two turbulence models, 
it is observed that the pressure difference in the SST k — o) 
model is slightly higher than the Standard k — £ model which 
is clearly observed in Fig. 8 also. 

3.2 Velocity Contours 

Fig. 11 and Fig. 12 depict velocity contours for the 
Standard k — e turbulence model and SST k — o) turbulence 
model, respectively, at different angles of attack. In these two 
figures, the effect of angle of attack is shown by changing the 
flow directions instead of rotating airfoil, which gives the same 
effect on airfoil. Similar velocity distribution pattern is 
observed for the Standard k — £ and SST k — co turbulence 
model, present in Fig. 11 and Fig. 12, respectively at lower 
angle of attack, whereas slightly varies in magnitude at higher 
angle of attack though patter of velocity distribution is similar. 
It is clearly observed from both figures that velocity above the 
upper is higher than the velocity below the lower surface. This 
happens due to the Bournollie’s principle. It is also observed 
from both figures that at lower angle of attack flow is attached 
to upper surface and lower surface as a result there is no zone 
of zero velocity which indicates flow separation happens. 
When the angle attack reaches to the value of 10® flow starts to 


separate, see Fig. 11, at the trailing edges which starts to move 
towards the leading edges at higher angles of attack. A similar 
pattern is also observed for the SST k — o) model, see Fig. 12, 
however, flow separation is identified comparative lower angle 
of at 8® compare to the Standard k — £ model. Due to this flow 
separation, lift coefficient starts to decrease after a certain angle 
of attack. 

4. Conclusion 

In this paper aerodynamic characteristics of NACA-4312 
airfoil have been studied using ANSYS Fluent by the Standard 
k — £ model and SST k — (o model turbulence model. At lower 
angles of attack both methods produce similar results however 
at higher angles of attack results vary in magnitude slightly 
though nature is similar. It is observed that the lift and drag 
coefficient increase with angle of attack though the lift 
coefficient shows a decreasing trend beyond the angle of attack 
13®. 
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